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Abstract 

' In this paper we describe the implementation of processes fifiZZ — > and f\f\HZ — > into the 

Oh, framework of SANC system. The /i stands for a massless fermion / whose mass is kept non-zero only 

in arguments of In functions and — > means that all 4-momenta flow inwards. The derived scalar 
' form factors can be used for any cross channel after an appropriate permutation of their arguments 

(s,t,u). We present the covariant and helicity amplitudes for these processes: for the former only 
in the annihilation channel /i/i — > ZZ, while for the latter in annihilation /i/i — > HZ and decay 
H — * Zfifi channels. We briefly describe additional precomputation modules which were not covered 
in the previous paper. For the processes /i/i — > HZ(ZZ) and decay H — ► Zfifi we present compact 
results of calculation of the accompanying bremsstrahlung and discuss exhaustive numerical results. 

As applications there are two types of the Monte Carlo generators for the process H — > 4/i. The 
first one is the generator based on a single resonance approximation for one of the Z bosons. The 
second one, exploiting the double resonance approximation, is not described in this article. For the 
generator in the single approximation we present a short description. 

Whenever possible, we compare our results with those existing in the literature. For example, 
we present a comparison of the results for H — > 4fj, decay with those obtained by MC generator 
Prophecy 4f. 
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PROGRAM SUMMARY 

• Title of program: SANC 

• Catalogue identifier. ADXK_vl_l 



• Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADXK_vl_l 



• Does the new version supersede the previous version?: Yes 

• Reasons for the new version: implementation of new processes; extension of an automatic generation 
of FORTRAN codes by the s2n.f package onto many more processes; bug fixes 

• Summary of revisions: 

— implementation of light- by- light scattering and Compton scattering with one virtual proton 
into QED branch 

— vast update of 2f2b node in EW branch 

— complete renovation of QCD branch 

• Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland 

• Designed for: platforms on which Java and FORM3 are available 

• Tested on: Intel-based PC's 

• Operating systems: Linux, Windows 

• Programming languages used: Java, FORM3, PERL, FORTRAN 

• Memory required to execute with typical data: 10 Mb 

• No. of bytes in distributed program, including test data, etc.: 

• No. of bits in a word: 32 

• No. of processors used: 1 on SANC server, 1 on SANC client 

• Distribution format: tar.gz 

• Nature of physical problem: Automatic calculation of pseudo- and realistic observables for various 
processes and decays in the Standard Model of Electroweak interactions, QCD and QED at one- 
loop precision level. Form factors and helicity amplitudes free of UV divergences are produced. For 
exclusion of IR singularities the soft photon emission is included. 

• Method of Solution: Numerical computation of analytical formulae of form factors and helicity 
amplitudes. For simulation of two fermion radiative decays of Standard Model bosons (W^, Z) and 
the Higgs boson a Monte Carlo technique is used. 

• Restrictions on the complexity: In the current version of SANC there are 3 and 4 particle processes 
and decays available at one-loop precision level. 



3 



• Typical Running time: The running time depends on the selected process. For instance, the symbolic 
calculation of form factors (with precomputed building blocks) for H — > e + e— Z process takes about 
10 sec, helicity amplitudes — about 10 sec, and bremsstrahlung — 1 min 10 sec. The relevant s2n 
runs take about 2 min 40 sec, 1 sec and 30 sec respectively. The numerical computation of decay 
rate for this process (production of bencmark case 3 Table) takes about 5 sec (CPU 3GHz IP4, 
RAM 512Mb, L2 1024 KB). 
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1 Introduction 



In this paper we continue to describe the computer system SANC Support of Analytic and Numerical 
Calculations for experiments at Colliders [T] intended for semi-automatic calculations of realistic and 
pseudo-observables for various processes of elementary particle interactions at the one-loop precision 
level. This is done in the spirit of the first description of the SANC system (see [I] and references therein 
which we recommend as a first acquaintance with the system). 

Here we consider the implementation of several processes of ffbb — > kind (where / stands for a 
fermion, b for a boson of the Standard Model (SM), while for concrete bosons we use A for the photon 
and Z, W^, H) One should emphasize also that the notation ffbb — > means that all external 4- 
momenta flow inwards; this is the standard SANC convention which allows to compute one-loop covariant 
amplitude (CA) and form factors (FF) only once and obtain it for a concrete channel by means of a 
crossing transformation. The present level of the system is realized in version v. 1.10. Compared to 
version v. 1.00, it is upgraded both physics-wise and computer- wise. As far as physics is concerned it 
contains an upgraded treatment of ud — > l + vi and du — > l~T>i processes (see Ref. [5]) and a complete 
implementation oi F — ► f + fi + f[ CC decays up to numbers and MC generators. (Here F and / stand for 
massive fermions and f\ and f[ for massless fermions of the first generation.) Although the version 1 . 10 
tree literally contains only t — > b+ 1 + + vi decay [3] , any decay of the kind F — > / + /l +/i may be treated 
in a similar manner and we are going to implement them into the next versions. The complete description 
of these CC decays will be given elsewhere [I]. Version 1.10 contains also the process H — > f\f\A in 
three cross channels [5] in EW branch, 77 — > 77 scattering [6] and 1 1 — * 77* in QED branch, as well as 
a new QCD branch [7]. 

New in version 1 . 10 are also several ffbb — -> processes, to whose implementation this paper is 
devoted. We describe here two of them: f\f\ZZ — > and f\f\HZ — > 0, the latter one being used in two 
channels — annihilation and decay. 

In the annihilation channel, these processes were considered in the literature extensively (see, for 
instance, [8] and [9]-[lQ]), however, we are not aware of publications devoted to the H — > decay. 

These processes are relevant for H search at LHC: the processes fxfi — * ZZ are one of the backgrounds 
while the one-loop calculations of the decay H — > Zfxfx was also used for an improved treatment of the 
decay H — > 4/j for an intermediate Higgs mass interval 130 GeV < M H < 150 GeV, see section [7] 

Furthermore, in the spirit of the adopted SANC approach, all 2f2b — > processes can be computed with 
off shell bosons thereby allowing their use also as building blocks for future studies of 5 — * processes. 

The process f\f\ZZ is very similar to the processes fflry (6 = 7, Z ,H) whose precomputation was 
described in detail in Ref. p] . Its tree level amplitude is represented by two diagrams in t and u channels, 



Fig. m 



fu=typeIU 



fu=typeIU 




i-channel 



u-channel 





fd=typeID 



a) fd=typeID b) 

Figure 1: Born f fHZ(ZZ) — > diagrams, t and u channels. 

5 



fd=typeID 
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Pi// H 



fu=typeIU 
Figure 2: Born ffHZ 



vd=typeFD 
diagrams, s channel. 
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Note, that in this and in the following figures, fu=typeIU etc. denote "types" of external particles, 
see Table 2 of [T] as well as the discussion in the beginning of section SJ 

For the process /i/i — ► HZ these two diagrams do not 
contribute in the tree approximation since fermion f\ is con- 
sidered to be massless. However, in this case there exists an 
s channel amplitude, Fig. O 

All ffbb processes are fully implemented at Level 1 of an- 
alytical calculations. Several new modules which compute 
the contribution of the bbb vertices to ffbb — > processes 
are added to the precomputation tree, as well as three other 
modules relevant for the s channel diagram. 

The modified "Precomputation" tree is shown in Fig. [5] 
and discussed in section [31 

The modified branch 2f2b for the "Processes" tree is 
shown in Fig. [3j It contains four new sub-menus /i/i — > 
ZZ, hh -> HZ, H -» hf x Z and H -> /ijU which in 
turn are branched into scalar Form Factors (FF) and He- 
licity Amplitudes (HA) (two for the process f\f\HZ corre- 
sponding to two annihilation and decay channels) and the 
accompanying bremsstrahlung contributions (BR). 

These processes are implemented also at Level 2, where 
the s2n.f package produces the results in the "Semi Ana- 
lytic" mode (see Fig. 20 of Ref. [I]). For the three decays 
we have relevant "Monte Carlo" generators which, however, 
are not yet implemented into the system. 
The paper is organized as follows: 

In section [2] we describe the covariant (CA) and helic- 
ity amplitudes for three of four new ffbb — > processes 
available in version 1 . 10. 

Section [3J contains a brief description of new precom- 
putation modules. In section [4] we describe in some more 
detail the renormalization procedure for the f\f\HZ — * 
process, i.e. calculation of FFs. Section [5] contains the results for the accompanying bremsstrahlung 
in the semi-analytic mode for two f\f\HZ — > channels. Section [S] contains numerical results for the 
processes /i/i — > HZ(ZZ) and decay H — > Zf\f\. Finally, section [7] contains a brief description a Monte 
Carlo generator for process H — > 4/i in the single resonance approximation. The first results of numerical 
comparison with those of Prophecy4f pT]-[T2] are also presented. 



Neutral Current 

□ f f -> A A 

□ A A -> f f 

□ f f -> H A 

□ f f -> Z A 
€S f1 f1 -> Z Z 

felt f -> Z Z (FF) 
fejf f -> Z Z (HA) 
v3 f1 f1 -> H Z 

Fwfl fl -> H Z (FF) 
Raf1 f1 -> H Z (HA) 
f1 f1 -> H Z (BR) 
Q H -> f1 f1 z 

fm H -> f1 f1 Z (FF) 
Rm H -> f1 f1 Z (HA) 
R?i H -> f1 f1 Z (BR) 

a h ->ti d a 

Fw H -> fl f1 A (FF) 
Fl23 H -> f1 f1 A (HA) 
□ Charged Current 



□ QCD 



Figure 3: New processes in the EW part. 
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2 Amplitude Basis, Scalar Form Factors, Helicity Amplitudes 



2.1 Introduction 

In this section we continue the presentation of formulae for the amplitudes of ffbb — > processes started 
in section 2 of Ref . [T] . As usual, we begin with the calculation of CAs corresponding to a result of the 
straightforward computation of all diagrams contributing to a given process at the onedoop level. It is 
represented in a certain basis of structures, made of strings of Dirac matrices and external momenta, 
contracted with polarization vectors of vector bosons. The amplitude is parameterized by a number of 
FFs, which we denote by T with an index labeling the corresponding structure. The number of FFs is 
by construction equal to the number of structures, however for the cases presented below some of the 
FFs can be equal, so the number of independent FFs may be less than the number of structures. For the 
existing tree level structures the corresponding FFs have the form 



T = 1 



a 



(1) 



where "1" is due to the Born level and T is due to the one-loop level. As usual, we use various coupling 
constants: 



e 
9 



My, 

M z 



etc. 



(2) 



Given a CA, SANC computes a set of HAs, denoted by H\ 1 \ 2 \ 3 ,.., where AiA 2 A 3 . . . denote the signs 
of particle spin projections onto a quantization axis. 



2.2 // -> ZZ process 

Here we present the CA of the process f{p2)J{pi) ~ > Z (pa) Z (pa) in the annihilation channel^, see Fig.[Tj 
It contains 10 left (7+) and 10 right (7-) structures: 



An 



k 



[i>3,i+{pi)^{pi)^^i "(M) +i>z"t + {pi) il {p 2 )uJ : i{s,t) 

+ fo7+(Pl)v(P2)nJ r 3~( S > t ) + fcl+{P2) v{P2)vFi {s,t) +^ 3 7 + ^ y J r 5 + (s,t) 
+ 7^37^7+-^6 + ( S ' l ) + 7m7+(Pi)^-^7 + ( S ) *) + l^l+(P2)u^s( S ^ *) 

+ 7v7+(Pi)^9 h ( s ! *) + 7»7+(P2) Al ^ r i l o(s, t]^u(p 2 ) ef,(p 3 )ef l (p 4 ) 
7+ -> Ft -> F- 



where 



Furthermore, 



and 



7± = I ± 75 • 



{pi +P2Y 



(P2+P3f = -t, {P2+Pif 



(3) 

(4) 
(5) 



lr The other channels are unphysical in this case. 
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Now we give the explicit form of the CA in the tree (Born) approximation: 



PfV (Pi) f 7^37^7+ + J lnl+(j>2)v + ^ ^ 3 7+<W 

-Jnl+(Pi)» +7^7+(pi)m +lvl+{P2)^j \u(p 2 )e^{p 3 )e^{p 4 ) 



'/ ~^ "S 



St, 7+ -> 7- 



(6) 



Note that this is decomposed into 12 structures of 20 and is highly asymmetric in t and u. This is due 
to our choice of the 4-momentum p 3 and of the ordering of Lorentz indices /x and v in Eq. © . 

Equation [5] may be parameterized by only two FFs if one introduces two "Born-like structures (BLS)" 
given by expressions in big round brackets by means of eliminating the 5th structure in favor of 

BLS; to this structure and to the corresponding (s, t) we assign the subindex "0": 



BLS,! - : - ( ~ + ~ ) 



-jl^l±{P2}u +7^7±(Pi)^ - 7^7±(Pi)p + 7^7±(P2%7) 



Moreover, between the 20 FFs there are four identities: 

Jf(a,t) = Jf{8,t), jf (s,t) = -Tf{s,t). 
Therefore, there are 16 independent FFs but 18 independent non-zero HAs for process /i/i — * ZZ: 



T2^ 7 +(s, t) ± ^ 8 +(s, t) T ^ 9 + (s, i) h 



fco^j^ Q + ~) J^ -(s,t) - |c±/3 [2JPf(s,t)- JP 2 -(s,t)- JP 3 -( S ,t)] 



+-±± 



T2^f (s, f) ± ^ (s, t) =p ^ (s, t) 

/3i /3: 



_ _+ 1 r - 

i u 



T+(s, t) - - sin 2 [2?+ (a, t) - T+{s, t) - F+(s, t)] 



+2 [(3 T F+(s, t) + cos d z Ff (s, t)] ± c T F+(s, t) ± c±F+{s, t)\, 



H. 



fc^ {281 



t u 



±0 



2 [/3 T ^" 6 -(s,i) +cos^f( S ,i)] T -c^M) T -CiJ^M)}, 



kf{2o 

s 



(¥L (i c + 2M 
t u s 







G- 


►=)] 



^ T /3 [T2cos^+( S ,t) T/£^ 2 + M) ±^ 3 +(s,t)] 
±/3 2 ^+( S , t) - ((3 C + T c T ) f+( a , t) T c T :F+(s, t) + /3i.F 9 + 0s, t)}, 



(8) 
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•0± 



kf\2a 



PI 0% 2M| (\ 1 

i M S \ t U 



H. 



+±o = 



- iC± /3 [±2cos^+(s,i) ±PLT+{s,t) (*,*)] 

M 2 ~\ 
±4-f-?+(8,t) - ((3% ± c±) - PLT£{8,t) t c±F+{s,t)\, 

PI p% 2M 2 (I 1' 



kf{25 2 el 
s 



t u 



t u 



T {s,t) 



+-c±p [T2co S i9 z ^(s,t)T(3-^2(^t)±P c + ^( Sl t) 



+o± 



TP\^ (s,t)- {P% ± c±) J-f (a, t) ± c±J^ (a, t) + /?^- (s, t) | , 

, + f r2 I"/?! /3i 2M 2 /l 1\1 , 
^ u s \t u ) J 

±^c T /?(2 cos «VT (*, t) + P1T 2 (s, t) - p%T^ (s, t)) 

M 2 1 



2 M 2 A T2 [ S 2 



P1_P% 

t u 









-=)] 



[2/^/£jf (a, t) + (/?L) 2 jf («, t) + (^) 2 jf(«, t)] 
+40-f-J$ (s,t)± 2p%T± (a, t) ± (s, t) - .F± (s, t)] j . 



Here we use the following shorthand notation: 

2 1 2 V2M ; 



p±l, 



Pl=P±cosd z , 



-c± , C± = 1 ± COSl? z , 



(9) 



(10) 



(11) 



and # z is the CMS angle between p 2 and p-$- The invariant i and the cosine cos$ z are related by 

t = M 2 Z - is(l-/3cosi? z ). 
The number 18 is the product of 2 initial massless helicity states and 3x3 states for the final Z bosons. 

2.3 /i/i -> #Z process 

There are six structures for the /i/i — > process if the fermion mass is neglected 



v(pi) ^l+^f^oi-^t) +p' 3 "f + (p 1 ) v J r f(s,t) + p , 3 - /+ (p 2 ) v J 7 £(s,t) )u(p 2 )e*(p 3 ) 



(12) 



The structures for the decay H — > f\f\Z may be obtained by simple replacement of 4-momenta 
— P3, P2 — > — P4, P4 — Pi (p3 — > P2) of the structures (fT2|) . 

Note, that the two terms ^^(s^t) correspond to the Born level. 

As far as HAs are concerned, we present them in both channels: annihilation and decay. 

2.3.1 HAs in annihilation channel f\f\ — > i?Z 



There are 6 HAs in 


this case: 






H+-+ 


k s c+ 


[fcf [f+(s,t)-F+(s,t)] -4a e ^ + (s,t)}, 




H-++ 


= — fcgC_ 


[fc+ [JT(«,t)-^(«,t)] +4<5 e ^ -( s ,i)}, 




H+— 


= — kyC- 


[fc+ "-^(M)] +4a e ^+( S ,i)}, 




H-+- 


k s c + 


[fcf [JT(«,t)- JT(«,t)] -46 e J^(s,t)}, 


F4cr e JP +( S ,t)}, 


T~t+-o 


k s k 2 < 




^-+0 


= ~k s k 2 < 


VA(s,M|,M2) [/^7/( s ,t) +/3iJF-( s ,t)] - 


f45 e JF -( S ,t)}. 



where 



k s = k 



f 



1 \/iM z 
s - MJ ' 



fc 2 = 



Mi 



\/2yfsM : 



M l ■ 

- sin ■&■> 



1 ± cos $ z , 
1 



/3i = /3±cosi? z , 



y/\(s,MlM*) 
8 + Ml -Ml ' 



Mi 



(a + Ml-M£)(l-0aoB0,). 



2.3.2 HAs in the decay channel H -> /1/1Z 

The six HAs in this case are somewhat different from the previous case: 



TL++- 
H-+- 



k s [k x [T+(s.t)-T+{s,t)} -4a f c + F+(s,t)}, 

k s {h [f+(s,t)-Ff(s,t)] + 4a f c-f+(s, t)], 
k s k 2 { y/\(M* , Ml sj [fi\T^ (a, t) + (i c _T 2 (s, *)] - 48 e F„ (s, t) } , 
{ - v/A(M2,M|, S ) (s, i) + /3L.F+ (s, t)] + Aa e T+ («,*)}• 



k s Q k 2 
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Here, 



fro 



1 ygMg 



h = y/\{M*,M*, a ) sin 2 i?/, 



(M 2 -Mj- S ) 



■ sin / , 



licosi?/, /3| = /3±costf/, /3 = 



VA(Mg,M|,.- 
ME - - ,■ 



(17) 



The number 6 is the product of 2 initial massless helicity states and 3 states of the final Z boson. 

Furthermore, s = is the invariant mass of the two fermions /, varying in the limits 4m 2 < s < 

(M H — M z ) 2 - and t is another independent kinematical variable, depending on s and an angle varying 
in the limits < #/ < ir 



t = Mi + 



Mi - Mi 



y/\(M*,Mi,s)cos0 } 



(18) 



The kinematical diagram of the process is shown in Fig. [4] 




Figure 4: H — > fxfxZ decay kinematics. 



The Higgs boson with momentum pi at rest, decays back-to-back into a Z boson with momentum p\ 
and a fermionic compound with 4-momentum p% + p^ and invariant mass s. This compound decays in 
its own rest frame into two back-to-back fermions with ■& t being the angle between p 3 in the compound 
rest frame and the direction of flight of the Z boson in the H boson rest frame. 
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3 Precomputation news 



The "Precomputation" tree of version 1 . 10 is shown in Fig.[5]with all modified sub-menus open, and all 
sub-menus closed which were not changed compared to version 1 . 00. In this section we briefly discuss 
what every new module does. 



□ 



Root 

9 SANC 
o- □ QED 

3 BN 

7 £3 Precomputation 
t- a Self 
t-'vS Boson 

¥ <3 Ren Blocks 

F* CalcBosRenConst 

CalcBosSEfunctions 
E?i CalcBosct 
ffHZ RenSelf 
F^i Boson Self 
Tadpole 
Q Fermion 
<3 Vertex 
a bbb 
*■ C3 Boson 
*- C3 Fermion 
<3 ffbb bbb 

ffAA bbb Vertex 
ffZA bbb Vertex 
FwffHA bbb Vertex 
ffZZ bbb Vertex 
ffHZ bbb Vertex 

<3 bff 

1m bff Vertex 
ffbb Vertex 
ffHZ bff Vertex 

*- □ Box 

> □ Processes 
* □ QCD 



Figure 5: New EW precomputation modules. 



First, we added a new folder accessible via menu 
sequence EW — > Precomputation — > Vertex — > 
bbb -> ffbb bbb with five modules fFXX bbb Ver- 
tex, XX=AA, ZA, HA, ZZ, HZ which compute 
three boson vertices of four topologies (see Fig. 11 of 
Ref. [TJ) to the corresponding ffbb processes, Fig. [6] 
The results of their calculations are saved to ffXX*.sav 
files to be loaded by corresponding modules computing 
FFs via chains EW — > Processes — * 4-legs — > 2f2b 
— * Neutral Current for the ffXX — > processes. 

These three boson diagrams contain both bosonic 
and fermionic components. The latter are precom- 
puted by the modules bbb Vertex in Boson and Fer- 
mion folders of the same level on the tree. Five mod- 
ules of ffbb bbb folder load them and then apply te- 
dious calculations involving in some cases the Schouten 
identity. There are many peculiarities in these calcu- 
lations, forcing us to have an individual module for 
each ffbb process. Note also that if the correspond- 
ing process has a Born-level s channel exchange as in 
Fig. [21 then the contribution of one-loop vertices is 
supplemented by the relevant counterterm cross. 

In the modules under discussion a summation over 
the exchanged boson B is performed. In general, four 
neutral bosons B = 7, Z, <f> and H can contribute if 
the fermion mass is not neglected, otherwise, only 7 
and Z contribute. 

For the processes ffAA -> 0, ffZA -> 0, ffHA -> 
the "left" bff vertex diagram, shown in Fig. [7] does 
not contribute, since in these cases the "right" vertex 
does not exist at the tree level. In general, and this 
is indeed the case for the processes ffZZ — ► and 
ffHZ — * 0, the "right" vertex exists at the tree level 
for B = Z,<fi°, therefore, the dressed "left" vertex has 
to be added to the precomputation tree. Note that it 
does not contribute for massless fermions if B = <j>° . 
For ffHZ — > this vertex is accessible via menu se- 
quence EW — > Precomputation — > Vertex — > bff 
-> ffHZ bff Vertex. For the ffZZ -> process, only 
B = H contributes, but then, for the massless /, the 
dressed "left" vertex vanishes. This is why we do not 
add the corresponding module in the latter case. 
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fd=typeID 



vu=typeFU 




Figure 6: "Right" bbb vertex in ffbb processes. 




Figure 7: "Left" bff vertex in ffbb processes. 



The presence of an s channel tree level diagram in the process f\f\H.Z — > (Fig. ^ forces us to take 
into account two more self energy diagrams, Fig. [5] 




fu=typcIU 



vu=typeFU fd=typcID 



fP2 

fu=typeIU 



Figure 8: Self energy ffHZ diagrams. 



vu=typeFU 
P*/H 




vd=typeFD 



The first one is accessible via menu chain: Self — > Boson — > Ren Blocks — > ffHZ Ren Self. Here 
Bi = 7, Z, cj) and B 2 = Z, again <fi°'s do not contribute for massless fermions. As will be explained 
in the next section, the second diagram is better to be combined with the "right" vertex, Fig. [S] 

Note that nothing is changed, compared to ffXA processes, as far as boxes are concerned. So, the 
Box sub-menu is as in version 1.00 pQ. 

The new modules contain calls to several new intrinsic procedures which will be described elsewhere. 
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4 Renormalizaton for ffHZ — > process 



In this section we describe how to use the FORM [T3] module which computes FFs for the process 
ffHZ — > 0. This description is supposed to help a user to understand the other modules computing FFs 
for any ffbb — > process. 

First of all, to use our basic declaration and notation we begin the file with 

#include Declar.h 
#call GlobalsO 

and define types of external particles, see Figs. [1] and O 



#ifdef 


'typelU' ; 


* 


<fu 


#ifdef 


'typelD' ; 


* 


f fd 


#ifdef 


'typeFU' ; 


* 


'vu 


#ifdef 


' typeFD ' ; 


* 


<vd 



#define typelDp "{2*( < typeID'°/ 2)-l+ < typeID'>" ; * 'fdp' 

Secondly, we fix four main steering flags to define the calculation Al scheme: 

1. define xi: xi = to test gauge invariance in Rg, or xi = 1 to work in £ = 1 gauge; 

2. define on: on = external photons are off mass-shell, or on = 1 photons are on mass-shell; 

3. define mf: mf = zero external fcrmion mass (i.e. pm('fd')=0), or mf = 1 it is not zeroed; 

4. define mp: mp = zero mass of the weak isospin partner of the fermion /, pm('fdp')=0, or mp = 1 
it is not zeroed. 

Actually, for the process under consideration, f\f\HZ — > 0, only the £ and mp definitions are mean- 
ingful since there are no external photons and we ignore the masses of external fermions throughout the 
calculations. But the mass of a weak isospin partner of an external fermion that appears in the internal 
loop may be kept nonzero. 

The ideology of building blocks (BB) is the key element for SANC development. The information about 
the main precomputed BB is stored in basic *.sav files (BSF). Note that for 4-particle ffbb processes 
all the BBs are 4-legs by construction. This trick will greatly simplify the procedure of projection of the 
covariant amplitude onto an independent basis of structures. 

Moreover, for the future development it is necessary to upgrade the database of the collected infor- 
mation area SANC: fields of program modules, fields of procedures and the bank of BSF. 

Any module computing FFs starts from loading of the calculated BB from the bank of BSF. These 
BSF contain the precomputed objects: self energies, vertices and boxes typically with off-shell bosons. 

They are precomputed not only to accelerate the calculations. Although all BSF may be, in principle, 
precomputed online, we remind that in some cases the CPU time for calculating off-shell boxes in 
gauge takes many hours, see section 3.4 of Ref. [lj. In such cases precomputation is strictly prohibited 
and the user must use already precomputed BSFs. 

We recall also that our precomputation procedure has indeed several levels; in the modules computing 
FFs we tend to use the results of the last level which contains already renormalized BBs: propagators 
and vertices, i.e. taking into account relevant counterterms and special vertices, [14J . However, they are 
full of residual UV poles and £ dependent terms, which cancel in the sum for a one-loop CA of a physical 
process. This is why we still use to word "renormalization" in connection with modules computing FFs 
rather than a simple "summation" . 

Typically, the loading of BSFs is organized in several steps. Let us consider the example of 
H -> fl fl Z (FF) module, see the tree in Fig. 
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• step self 

Here we manipulate objects from BSF ffHZSelfschxi'xi"fu"fd"vu"vd'.sav 

We extract from its volume the BB of bosonic self energy in the s-channel, BSEsch'fu"fd"vu"vd' see 
left diagram in Fig. [5J 

• step vertex 

Moving further over the renormalization procedure at this step we load three BSFs: 

ffHZVertbffxi'xi"fu"fd"vu"vd'.sav; 

ffHZVertbbbxi'xi"fu"fd"vu"vd'.sav; 

ffbbVertxi'xi'on'on'mf'mf 'mp'mp"fu"fd"vu"vd'.sav. 

From these BSFs we extract various types of vertices correspondingly: 

— VertBff'i"fu"fd"vu"vd' with i=l,2,3,4 standing for £ A , £ z , £ w , and no £ vertex clusters originating 
from the diagram of Fig. [2 

— Vertbbbbos'fu"fd"vu"vd' and Vertbbbfer'fu"fd"vu"vd' — the bosonic and fermionic components 
of three-boson vertices shown in the diagram Fig. [6j where the former contains counterterms, 
the special vertex and the right diagram of Fig. [8j These tadpoles cancel the £ z dependence of the 
three-boson vertices, giving an opportunity to assign this contribution to i=3. Finally, the fermionic 
component should be naturally assigned to i=4; 

— abelian VertT'i' and non-abclian vertT'i' vertex clusters in t and u channels l=t,u with cluster 
index i=l,2,3,4 and k=22,33,24,42,44, see section 3.4.2 of Ref. Q] for a description of the latter. 

• step boxes 

Here the most complex building blocks — off-shell boxes are loaded from four BSFs: 
ffbb3Tlxi'xi'on'on'mf'mf 'mp'mp"fu"fd"vd"vu'.sav; 
ffbb3T3xi'xi'on'on'mf'mf 'mp'mp"fu"fd"vu"vd'.sav; 
ffbb22T5xi'xi'on'on'mf'mf 'mp'mp"vd"fd"vu"fu'.sav; 
ffbb33T5xi'xi'on'on'mf mf 'mp'mp"vd"fd"vu"fu'.sav. 

They contain precomputed boxes of topology Tl from the expression S3Tl'xi"on"mf "mp"fu"fd' 
'vd"vu', the boxes of topology T3 from the expression S3T3'xi"on"mf "mp"fu"fd"vu"vd', and of 
topology T5 with cluster index kl = 2,3 , i.e. with virtual Z and W bosons from the expression 
S'kl"kl'T5'xi"on"mf "mp"vd"fd"vu"fu', see section 3.4.2 of Ref. pQ. 

Only those box topologies and clusters are loaded which give a non-zero contribution for m/ = 0. 

• step Sum 

Finally, we sum all contributions. Four expressions, Sum'i', corresponding to cluster index i=l,2,3,4 
are being constructed here. The first three of them may carry only one gauge parameter each, the 
latter carries none. 

After construction of four Sum'i', the module continues with various kinds of transformations (in 
particular, involving an algebra of Gram determinants) which prove the cancellation of gauge parameter 
dependences in first three Suml,2,3 and the cancellation of the residual UV poles between FFs with 
cluster indices i=3 and i=4. 

After the comment Preparing Structures for obtaining FFs: the six basis elements of the CA, 
Eq (fT2|) . are created and the 6x4 FFs are projected out of this CA. 

The final step is formatting of the BSF FFflflHZ.sav with 24 FFs — FFgp'k'T and FFgm'k'T 
(k = 0, 1, 2, i = 1, 4) — for subsequent processing by s2n.f software. 
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5 Bremsstrahlung in fifiHZ — > processes 

In this section we present the list of short final results for the contribution of accompanying bremsstrahlung 
processes. 



5.1 Bremsstrahlung in fifi — > HZ annihilation channel 

The tree level diagram of this channel is shown in Fig. [2] The corresponding total Born cross section 
reads: 



Ror 



Gp(vj + a 2 ) M t\J H s > M z> M h) 



12tt (M| - s) 2 + M|r| 
There are only two initial state (ISR) bremsstrahlung diagrams: 



\ + \ ~ M l) + ^2 « - M D 



(19) 





Figure 9: Bremsstrahlung diagrams in annihilation channel 
QED corrections due to virtual and soft photons are proportional to the Born cross section: 
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(20) 
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7T 



In 
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In 
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(21) 



with infrared divergence In A 2 being canceled ou10. The hard photon contribution of the differential cross 
section in s' = —(p3 + pa) 2 has the following factorization property: 



da 



Hard 



ds' 



ir * s 2 (s — s') 



,s 2 + s' 2 



111 



C7 Born (s') 



(22) 



It may be integrated over s' leading to a rather compact expression for cr tot (see relevant module at SANC 
tree), yielding a lnw which cancels against corresponding term in Eq. (|21l) . 
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h 

Figure 10: The H — ► f\f\Z decay tree level diagram. 



5.2 Bremsstrahlung in H — > f\f\Z decay channel 

Here we consider the decay channel H — > j\f\Z. We begin with the tree level diagram, Fig. [TOl 

The corresponding tree level double differential width, depending on two kinematical variables s, ??/ 
discussed in section 12.3.21 and with kinematics shown in Fig. SI reads 




Figure 11: The H — > fifiZ decay, bremsstrahlung. 



2 Note, that the "Soft" contribution to the process fifi — * ZZ is also described by Eq J21l l. 
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The fully differential phase space is characterized by hve kinematical variables which we choose as 
follows: 



Z7T Z7T 



(24) 



where r = — {pa + P5) 2 is the lepton-photon invariant mass. 

The 3-step kinematical cascade develops as a sequence of three 2-body decays shown in Fig. [T2l 




Figure 12: Bremsstrahlung kinematics, 
with three corresponding two body phase spaces: 
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(25) 



The gauge invariant QED part of the complete one-loop EW correction is subdivided into virtual, soft 
and hard photon contributions. The virtual one comes from the two Born-like FFs with cluster index i=l. 
It is proportional to the Born width and contains the infrared divergence parameterized by the photon 
mass A: 



2-pVirt 



d 2 r 



j2t^Boiii 



ds d COS $ 



f 




ds d cos 1? 




1 



In 



3 -In 



-4Li 2 (l)-2 



(26) 

The soft photon contribution is also proportional to the Born one; its infrared divergence cancels against 
the virtual contribution. It contains also a logarithm with soft-hard separator u>: 



■2 r Soft 



d 2 r 



'2-pBorn 



d 2 r 



ds d cos i9f dsd cos ■df n 




In 



A 2 



2 In 



2Q 
m f 



In 



In 



1 



Li 2 (1) + 1 



: . . (27) 

The hard photon contribution after integration over three kinematical variables dtp 4 , d cos §4 and dr (the 

first two vary together in full angular Ait limits and m 2 < t < (a/s — m f) 2 ) is 
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'2^Hard 



d 2 r 



2-nBo 



d 2 r 



ds d cos ■d / 



ds d cos 7T ^ 



-3Lia(l 



-Ql -2 In — 




(28) 



The total QED correction, sum of above three contributions, is free not only of infrared divergence 
and of soft-hard separator, but also free of final fermion mass singularity in accordance with the KLN 
theorem [TSj-fTo]. 



2r^Total 



d 2 r 



2t^Bo 



ds d cos 



f 



rl 2 Y 

— -Q% 

ds d cos df it 1 



k B Q 
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Finally, if one integrates over dcost?/, the well known Z decay correction factor restores: 



dY 



Total 



ds ds 

Therefore, the QED part of the correction is small, 



1 

0.2% 



4 7T J 



(29) 



(30) 



6 Numerical results and comparison 

In the numerical calculations by s2n package we use two precompiled libraries: SancLib_vl . 00 and 
looptools 2.1 [17]. 



6.1 Numerical results for Electroweak corrections 
6.1.1 Process f x fi ~> HZ 

For this process we present in Table 1 the results of a tuned triple comparison of the one-loop electroweak 
corrections, excluding the gauge-invariant QED subset of diagrams (vertex and electron self-energy) and 
real bremsstrahlung. The input parameters are taken as in |18j . Table 1 shows 6-7 digits agreement 
between the three calculations. 



v/s, GeV 


Mh, GeV 


.18J 


Grace-Loop 


SANC 


500 


100 


4.1524 


4.15239 


4.15239 


500 


300 


6.9017 


6.90166 


6.90166 


1000 


100 


- 2.1656 


- 2.16561 


- 2.16560 


1000 


300 


- 2.4995 


- 2.49949 


- 2.49949 


1000 


800 


26.1094 


26.10942 


26.10942 


2000 


100 


-11.5414 


-11.54131 


-11.54136 


2000 


300 


-12.8226 


-12.82256 


-12.82256 


2000 


800 


11.2468 


11.24680 


11.24680 



Table 1: Comparison of percentage correction to the total cross section e + e — * ZH between Ref. [18j . 
Grace-Loop [llj and SANC. 

Beside the input given in [18], Mw is crucial for a precise comparison. The following Mw masses 
have been used: M w = 80.231815GeV (M H = lOOGeV), M w = 80.159313 GeV {M H = 300GeV), 
Mw = 80.081409GeV (Mh — 800GeV), following A. Dcnner private communication, as referred to 
in [TO]. 
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6.1.2 Process /i/i -> ZZ 

For this process we compared only "virtual + soft" corrections in the conditions of Tables 1,2 of Ref.[5] 
with input parameters tuned carefully. In this case we do not find a good agreement: SANC numbers 
happened to lie about 10% lower. We shall go back to searching for the origin of this discrepancy after a 
tuned comparison of the very similar process /i/i — > ZA with the results of Ref. [217] . The implementation 
of the latter process into the SANC system is nearly finished. 

6.2 Numerical results for real and complete corrections 
6.2.1 Hard bremsstrahlung in /i/i — > HZ annihilation channel 

In Table [5] we present typical results of a triple comparison of the Born cross section and the cross section 
of hard photon bremsstrahlung between two calculations within SANC (semi-analytic, Eq|22l and MC) 
and those of CompHEP for E-y > lGeV. Here we used m_ff=130 GeV and the other parameters as in 
CompHEP. 





a, pb 


Vs, GeV 


250 


300 


500 


1000 


2000 


Born (SANC) 
Born (CompHEP) 


0.21984(1) 
0.21984(1) 


0.17454(1) 
0.17454(1) 


0.056890(1) 
0.056889(1) 


0.012898(1) 
0.012898(1) 


0.0031322(1) 
0.0031322(1) 


Hard (SANC, s2n) 
Hard (SANC, MC) 
Hard (CompHEP) 


0.080309(1) 
0.080307(1) 
0.080306(2) 


0.091168(1) 
0.091166(1) 
0.091168(2) 


0.043650(1) 
0.043649(1) 
0.043651(1) 


0.013246(1) 
0.013246(1) 
0.013242(3) 


0.0040293(1) 
0.0040293(1) 
0.0040287(4) 



Table 2: Comparison of the Born cross section and hard photon cross section of e + e — > HZj reaction 
for E-y > 1 GeV. 

The two SANC results perfectly agree within statistical errors. The agreement with CompHEP also 
looks quite good. 



6.2.2 Hard bremsstrahlung in fifi — » ZZ annihilation channel 

In Table [3] we present the results of a double comparison of the Born cross section and the cross section 
of hard photon bremsstrahlung between the calculation within SANC (MC) and those of CompHEP for 
E-y > lGeV. Here we used run =130 GeV and the other parameters as in CompHEP. 



a, pb 


Vi,GeV 


250 


300 


500 


1000 


2000 


Born (SANC) 
Born (CompHEP) 


1.0758(1) 
1.0758(1) 


0.82971(1) 
0.82971(1) 


0.40644(1) 
0.40644(1) 


0.14815(1) 
0.14814(1) 


0.049760(1) 
0.049761(1) 


Hard (SANC, MC) 
Hard (CompHEP) 


0.52262(1) 
0.52258(2) 


0.46800(2) 
0.46801(1) 


0.29114(1) 
0.29111(2) 


0.13375(1) 
0.13374(1) 


0.54563(1) 
0.54563(6) 



Table 3: Comparison of the Born cross section and hard photon cross section e + e —* ZZj reaction for 
Ey > 1 GeV. 

Again, we see a very good within larger statistical errors in CompHEP. 

3 For this channel we don't have the semi-analytic result. 
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6.2.3 H -> fifiZ decay channel 

For this decay we present the complete one-loop correction. 

We present numbers, collected in the scheme for the standard SANC INPUT: PDG(2006) [21] 



G F 


= 1.16637- 10- 5 GeV~ 2 , 


a(0) 


= 1/137.03599911, 


a s (M z ) = 


0.1187, 


M w 


= 80.403 GeV, 




= 2.141 GeV, 






M z 


= 91.1876 GeV, 


r z 


= 2.4952 GeV, 






M H 


= 120 GeV, 










m e 


= 0.51099892- 10- 3 GeV, 


m u 


= 62 MeV 


m d = 


83 MeV, 




= 0.105658369 GeV, 


m c 


= 1.5 GeV, 


m s = 


215 MeV, 


m T 


= 1.77699 GeV, 


771ft 


= 4.7 GeV, 


m t = 


174.2 GeV, 


T t 


= 1.551 GeV. 











The only exception is the Higgs boson mass for which we again take M H = 130 GeV in this section. 
Table 0] shows the double and singly differential decay width for the decays H — ► e + e~Z and H — > 
for a set of s and cost?;. 

From Table [4] it is seen that at the edges of cos di and near the fermionic threshold the double 
differential width shows a 1/s behaviour, typical for Coulomb interaction. The origin of the Coulomb 
peak at the one-loop level may be easily understood. First we note that H — * width does not vanish 
for an on-shell photon with Q 2 = 0, see first Fig. [121 




h 

Figure 13: H — > Z7 decay and Compton singularity. 

Therefore, the one-loop amplitude for H — > Zf\f\ with virtual photon exchange will show a ~ 1/s 
behaviour (with s = —Q~ ( )- This, in turn, will lead to the ~ 1/s behaviour of both the double and single 
decay differential widths. This conclusion is fully confirmed by the numbers in Table [H 

Recalling now the limits of s, 4m 2 < s < (M H — Mz) 2 , one might expect the appearance after 
integration over s of the big logarithm ln((M H — M z ) 2 /mj), with a final state fermion mass singularity 
However, the 1/s region is very narrow and it is largely washed out not only by a soft cut on the variable 
s but even by the plain integration over s. 

Finally let us discuss the total width for muon channel. In Born approximation it is r Born = 5.592- 10 -6 
GeV, while with the complete EW corrections it is r Born+1 - loo P = 5.774 • 10~ 6 GeV. So, the correction 
in scheme amounts to 3.2%. 
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H -» e+erZ 


Part 1, d?T/ds d cos tf; • 10 s , GcV" 1 




y/s, GeV 


1 


3 


10 


38 


Born 






0.04760 


0.04910 


0.06596 


0.09604 


1-loop 


COS^; = 


±0.9 


0.54806 


0.10413 


0.07078 


0.09875 








10.5123 


1.12053 


0.07301 


0.02829 


Born 






0.18725 


0.18788 


0.19473 


0.09768 


1-loop 


COS^; = 


±0.5 


0.53457 


0.22763 


0.20006 


0.10044 








1.85478 


0.21159 


0.02737 


0.02823 


Born 






0.24960 


0.24983 


0.25221 


0.09842 


1-loop 


COsfli = 


0.0 


0.52882 


0.28303 


0.25801 


0.10120 


5 






1.11867 


0.13289 


0.02296 


0.02820 


Part 2, cZT/ds • 10 y , GcY" 1 




Born 




3.330 


3.345 


3.511 


1.949 




1-loop 




10.73 


4.186 


3.618 


2.004 


1-loop/Born 




2.224 


0.252 


0.030 


0.028 




Part 1, d^T/dsdcostfi ■ 10 s , GeY -1 


Born 






0.05533 


0.04996 


0.06602 


0.09603 


1-loop 


COS??; = 


-0.9 


0.54437 


0.10480 


0.07081 


0.09875 


S 






8.83770 


1.09736 


0.07265 


0.02829 


Born 






0.18573 


0.18770 


0.19470 


0.09768 


1-loop 


COS^; = 


-0.5 


0.52536 


0.22739 


0.20003 


0.10044 


5 






1.82850 


0.21142 


0.02740 


0.02823 


Born 






0.24395 


0.24919 


0.25214 


0.09842 


1-loop 


COsfii = 


0.0 


0.51713 


0.28238 


0.25795 


0.10044 


S 






1.11982 


0.13317 


0.02303 


0.02823 


Born 






0.18573 


0.18770 


0.19470 


0.09768 


1-loop 


COS??; = 


0.5 


0.52538 


0.22739 


0.20003 


0.10044 


5 






1.82862 


0.21144 


0.02740 


0.02823 


Born 






0.05533 


0.04996 


0.06602 


0.09603 


1-loop 


COS$i = 


0.9 


0.54441 


0.10480 


0.07081 


0.09875 


S 






8.84841 


1.09749 


0.07266 


0.02829 


Part 2, dT/ds • 10 9 , GcY" 1 




Born 




3.327 


3.344 


3.511 


1.949 




1-loop 




10.57 


4.184 


3.617 


2.004 


1-loop/Born 




2.176 


0.251 


0.030 


0.028 



Table 4: The double and single differential widths for two decay channels: H — > e+e~ Z and H — ► fi + n~ Z. 
The table contains two parts for each channel. Part 1: first row: the double differential decay width 
d 2 r /ds d cos $i ■ 10 8 GcV _1 at the Born level; second row: the double differential decay width at the 
1-loop level; third row: relative correction S = d 2 r 1_loop /c? 2 r Born . Numerical values are truncated to 6 
figures. Part 2: the same set for the single decay differential width dT/ds - 10 9 , GeV -1 . Numerical values 
are rounded. 
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6.2.4 Hard photon radiation in H — ► fifiZ decay 

Here we present the triple comparison of the hard photon bremsstrahlung contribution with a cut on 
the lepton-photon invariant mass r = — (j?4 + p$) 2 again between two calculations of SANC and one of 
CompHEP. 



T • 10'', GeV 


yj Train 




1 


2 


3 


5 


10 


e 


s2n 


57.456 


17.982 


13.094 


10.395 


7.2097 


3.4621 


MC 


57.5(2) 


17.981(3) 


13.098(2) 


10.398(1) 


7.2100(5) 


3.4619(7) 


CompHEP 


55.3(5), unstable 


17.96(1) 


13.10(1) 


10.39(1) 


7.198(5) 


3.458(2) 


fi 


s2n 


18.917 


9.6862 


6.8385 


5.3224 


3.5885 


1.6396 


MC 


18.916(3) 


9.6869(6) 


6.8384(5) 


5.3227(4) 


3.5886(8) 


1.6399(3) 


CompHEP 


18.91(1) 


9.678(7) 


6.837(5) 


5.323(4) 


3.586(3) 


1.646(1) 


T 


s2n 


5.9714 


out of 


5.0085 


2.9539 


1.7377 


0.68272 


MC 


5.9717(4) 


kinematical 


5.0081(5) 


2.9542(3) 


1.7375(1) 


0.68269(5) 


CompHEP 


5.983(4) 


region 


5.016(4) 


2.959(2) 


1.741(1) 


0.6832(5) 



Table 5: The decay width T in GeV, massive case, as a function of ^r rnin from SANC s2n, SANC MC and 
CompHEP 

In Table[5] 'V Tmin |is=o jGeV" denotes T m i n derived by the formula r m i„ = mf + 2miE 7iTn i n . 
As seen, two SANC numbers agree within MC errors and there is a reasonable agreement with CompHEP 
everywhere but upper left corner (soft radiation by electrons) where CompHEP show a tendency to be 
unstable. 



7 A Monte - Carlo generator for H — > 4/i 



A Monte Carlo generator of unweighted events for process H 
application of the building blocks ideology of SANC. 
In our generator we implement two building blocks 
at the one-loop level: H — > f fZ(pf) and Z — > 
ff(l)- We merge these two blocks and create a 
link between them by means of the Z* (resonating 
boson) line with Brcit-Wigner mass distribution. 
In this spirit we create the generator in the single 
resonance approximation. We also built a double 
resonance generator, where we incorporated reso- 
nance approximation in two Z* lines. 

The range for application of the single reso- 
nance approximation was found to be 120 GeV< 
M H < 160 GeV and for the double resonance ap- 
proximation M H > 180 GeV. 



4/x is the first example of a possible 




Figure 14: H — » 
approximation. 



4/i in the single resonance 
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These conclusions are illustrated by Fig. [15] where we present the results of calculations of the tree 
level width of the decay H — ► 4/i for the three cases: 

1) solid line: results of complete tree level calculations neglecting effects of identical final state muons; 

2) dash dotted line: single resonance approximation, 



r 



H->2fj,Z^Z^2fi, 



3) dashed line: double resonance approximation, 



^2— res 



r 2 

1 7, 



(31) 



(32) 



le-02 



le-03 - 



le-04 - 



le-05 - 



le-06 - 



le-07 = 



le-08 - 



le-09 




100 



150 



200 



M H , GeV 



250 



Figure 15: Born width of the decay H — » 4/z in three approximations. 



The loop-corrected result is the linear combination of three types of events: Born with Identity events 
minus Resonanse Born events plus One-loop Resonance events. We feel that this is the main shortcoming 
of the generator. 

Let us consider each type of events: 

• "Born with Identity" events means a branch which computes the distributions without radiative 
corrections but with effects of identical muons; 

• "Resonance Born" events means resonance approximation for one of the Z bosons, i.e. H — > ZZ* — > 

where Z* is the resonating Z. Here, the two building blocks in Fig. [14] are calculated 
at the Born level; 

• "Resonance One-loop" events is implemented in the same spirit as "Resonance Born" , but building 
blocks H — > ffZ(-f) and Z — > //(7) are calculated at the one-loop level. 



24 



The codes of both MC generators can be obtained from the authors by request. 

Recently appeared a new MC code Prophecy4f, see Refs.|llj-|12j. realizing calculation of the complete 
one-loop corrected partial widths of the H — > Al channels. We present a preliminary comparison between 
MC Prophecy 4f and SANC in Table H 



y/s, GeV 


120 


130 


140 


150 


160 


Prophecy4f 


7.053(3) • 10" 8 


2.3769(9) • 10- y 


6.692(2) • 10- y 


1.6807(6) • 10- B 


4.006(1) ■ 10- B 


SANC (G„) 


7.197(3) • 10-" 


2.4079(8) • 10~ Y 


6.743(2) • 10-' 


1.6842(5) • 10- b 


3.962(2) • 10- b 


S,% 


2.04 


1.01 


0.76 


0.21 


-1.10 


SANC {a) 


6.938(2) • 10" 8 


2.343(1) • 10- y 


6.594(2) • 10- y 


1.6534(5) • 10- b 


3.915(1) • 10- b 



Table 6: Comparison for partial width for decay H — » 4/i in G M scheme for M H = 140 GeV between 
Prophecy 4f and SANC. 



As seen from the Table, there is ±1% agreement in the mass range 130-140 GeV, degrading at the 
edges of the interval [120-160], that finds its natural explanation in Fig[15l Moreover, Prophecy4f uses 
the complex-mass scheme and takes into account several higher order corrections. One has to emphasize, 
however, that SANC calculations in a and G M schemes differ by about 2%. This can be considered as 
a rough estimate of the theoretical error. Prophecy4f numbers lie basically inside the range of SANC 
predictions. 

The generator in the single resonance approximation described in this section was used for a MC 
simulation of H — ► 4/i decay in the ATLAS detector and the results were compared with simulation by 
PYTHIA, showing notable deviations, see [52]. This fact demonstrates the importance of higher order 
corrections and the necessity to reduce the theoretical error. 

8 User Guide 

8.1 Benchmark case 3: the process H — ► f\f\Z 

Here we consider the 2/26 NC process H — > f\f\Z. 

One can open the relevant branch of the SANC tree as follows: 

EW -> Processes -> 4 legs -> 2f2b -> Neutral Current -> H^flflZ 

For this process there are three FORM programs: (FF) Form Factors, (HA) Helicity Amplitudes, and 
(BR) Bremsstrahlung. Each of them in turn is opened, compiled and run as described in Section 6 of 
Ref.pQ. 

For the process H e + e~Z we have in the Console window the particle indices shown in Table[7J 
Table 7: Assignment of particle indices for the process H — > e + e~ Z 



typelU = 


4 


initial partile (H-boson) 


typelD = 


12 


final particle (electron) 


typeFU = 


12 


final antiparticle (positron) 


typeFD = 


2 


final particle (Z-boson) 
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These can be changed to typelD (typeFU) = 13,14 for up- and down-quarks in the hnal state of the 
processes H — > (uu, dd)Z by editing the particle numbers as explained in Section 6 of Ref. [T]. @. 

Next bring the Fortran Editor sheet of the Editors List and the Numeric Form panel to the 
foreground. Shown in the Numeric Parameter sheet are the particle masses in GeV and the invariant 
mass of ffZ compaund in GeV, also the cosine of the angle $i defined in FigHJ 

Click on the Rehash button at the bottom of the Numeric Form panel: the main module of 
FORTRAN code appears in the Fortran Editor sheet of the Editors List. Then click on Compile. The 
final answer appears in the Output field. It consists of the parameters used (a, Gf, particle masses, the 
't Hooft scale /i and the invariant mass of compaund), and the resulting differential width d 2 T/ds d cos i?j 
in the Born approximation and Born+one-loop. The results for the default parameters and for several 
scattering angles are summarised in Table [SJ 

Table 8: The double differential widths for H — > e + e~Z channel in a-scheme: first row: the double 
differential decay width d 2 T / ds dcosdi ■ 10 8 GeV~ 1 at the Born level; second row: the double differential 
decay width at the 1-loop level; third row: relative correction 8 = <i 2 r 1 ~ loop /(i 2 r Born . Numerical values 
are truncated to 6 figures. 



H -» e + e~Z 


Part 1, d^T/dsdcostii ■ 10 s , GeV -1 




v/i, GeV 


1 


3 


10 


28 


Born 




0.02019 


0.02144 


0.03505 


0.04261 


1-loop 


costf; = ±0.9 


0.21060 


0.04321 


0.03874 


0.04602 


5 




9.43022 


1.01508 


0.10537 


0.07984 


Born 




0.07914 


0.07964 


0.08478 


0.04353 


1-loop 


costf; = ±0.5 


0.21495 


0.09898 


0.09150 


0.04701 


S 




1.71589 


0.24281 


0.07922 


0.07976 


Born 




0.10546 


0.10562 


0.10698 


0.04394 


1-loop 


costf ; = 0.0 


0.21695 


0.12394 


0.11510 


0.04745 


S 




1.05716 


0.17343 


0.07586 


0.07972 



Input parameters can be changed by editing the appropriate field of the Numeric Form panel and 
pressing the Rehash button. Again the Rehash button must be pressed before pressing Compile. 

To produce whole Table [8] one can set flag tbprint = 1 in the Fortran Editors sheet. After editing 
the code just press Compile; there is no need to press the Rehash button. 

One can also produce the differential width dT/ds and total width T in GeV by integrating the above 
differential width. To produce these numbers one can set flag inf lag = 1,2, respectevely, in the Fortran 
Editors sheet. After editing the code just press Compile; again, there is no need to press the Rehash 
button. 
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